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ABSTRACT 

Because underwater gravity surveys are usually conducted in coastal waters, 
several unique corrections are required to reduce the observed data for calculation 
of free-air and Bouguer anomalies. The most logical sequence for the application 
of corrections would be: terrain correction (if necessary) to reduce bathymetric 
irrequiarities to a flat, horizontal bottom; removal of the upward attraction of 
the water above the meter; free-air correction to bring the meter measurement to 
the surface of the reference spheroid; filling in the ocean below the spheroid with 
rock for a Bouguer anomaly calculation or with water for a free-air anomaly 
calculation. The terrain correction must be handled carefully as the bathymetric 
irregularities are immersed in water of finite density and, due to the close 
proximity of land, topographic irregularities immersed in air may be within the 
range of the correction graticule grid. Graticule compartment weighting factors 
for making this correction may be readily determined by breaking the bottom- 
land surface profile into several segments delimited by the gravimeter depth and 
the sea surface . 

This work was supported by: Office of Naval Research through the Naval 
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CORRECTIONS FOR UNDERWATER GRAVIMETRY 



INTRODUCTION 

The concept that present shorelines represent generally insignificant 
discontinuities to the stratigraphy and structure of crustal rocks found on 
land and the adjacent continental terrace is basic to geology. However, the 
problems of safe ship navigation, sampling outcrops and crustal rocks buried 
by sediments, and conducting continuous geophysical surveys in the shallow 
coastal waters present difficulties in the testing of this continuity. For 
example, the correlation of land gravity measurements with sea-surface gravity 
surveys often requires interpolation of data across 5 to 10 km of unsurveyed 
terrain between the shoreline and shiptracks. One technique for filling this 
void is through the use of the underwater (bottom) gravimeter from a shallow- 
draft vessel. Although underwater gravimetry requires laborious field work, 
the fact that the meter is placed on the ocean bottom increases the accuracy 
of the observed values over those measured at the sea surface by avoiding the 
problems of removing the ship motion. Also, in coastal waters sea-surface 
gravity data may be unusable due to the degredation caused by severe ship mo- 
tion and excessive course changes to avoid navigational hazards. 

Although most discussions of underwater gravity surveys include some men- 
tion of data corrections, there is a general dearth of detail for a person con 
fronted by this situation for the first time. In addition, the terrain correc 
tion for bottom gravity data can become an important factor along continental 
shelves cut by steep submarine canyons and valleys; the application of this 
correction to data collected on the ocean bottom raises some important pro- 
blems to the investigator. The object of this paper is to present a simpli- 
fied discussion of data reduction applicable to underwater surveys. Tradi- 
tionally, gravity corrections are discussed in the following order: free-air, 

Bouguer, and terrain. This order is significant in terms of increasing diffi- 
culty of calculation, but it is not necessarily the most logical order for 
understanding the physical significance of these corrections and their proper 
application to a specific geologic region. 

Although underwater gravity surveys have been conducted from a diving 
bell, lowered to the ocean bottom with a gravimeter and its operator inside 
(Frowe, 1947), the most practical method employs a gravimeter encased in a 
pressure housing and controlled remotely from the ocean surface through elec- 
trical conductors in the lowering cable. Pepper (1941) and LaCoste (1967) 
discussed the functioning of such gravimeters and Beyer et al. (1966) deter- 
mined the feasibility and precision of measurements made with a modified 
LaCoste and Romberg gravimeter in water depths up to about 900 m. At the 
Naval Postgraduate School, the LaCoste and Romberg Model H6G Underwater Gravi- 
meter is used for shallow-water surveys. This instrument is essentially a 
land meter mounted inside a shell of two thick aluminum hemispheres supported 
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by a triangular base. Regulation of the temperature of the meter is accom- 
plished with a thermister-transistor circuit in the meter housing within the 
pressure shell. Clamping and unclamping of the gravimeter beam or weight, 
leveling of the meter within the shell, and manipulation of the measuring 
screw are all remotely controlled from the sea surface. The maximum bottom 
slope permissible for remote leveling of the meter is 15°. From the ship, the 
meter is lowered and raised by a hydraulic winch containing a slip ring assembly 
for electrical continuity between the meter and the control box. The control 
box contains functions for meter temperature control, leveling, beam clamping 
and unclamping, and measuring screw control. Also on the control box are indi- 
cators showing possible flooding of water into the meter, unacceptable tilting 
of the meter due to excessive bottom slope or tipping of the meter on the bottom, 
and the water depth as measured by a pressure sensor on the meter, calibrated 
for salt water. In 50 m of water, the time required for a measurement is approxi- 
mately 4 minutes. The accuracy of the meter, determined by periodic checks at 
a base station, is between 0.1 and 0.2 mgal (La Coste, 1967; Brooks, 1973), com- 
pared with accuracies of 1 to 2 mgal possible for sea surface gravimeters when 
good navigational control is available and a reasonably calm sea state prevails 
(Wall, Taiwan!, and Worzel, 1966; Beyer et al. , 1966: Dehlinger and Chiburis, 
1972). 

CORRECTION MODELS 

To demonstrate the various corrections to be applied to the observed under- 
water gravity values, the following models will be described: flat, horizontal 

bottom: smooth sloping bottom; irregular bottom. General equations for the 
corrections are given along with specific formulas assuming densities for sea 
water and crustal rocks. These formulas may be modified for particular loca- 
tions using more appropriate densities. 



Flat, horizontal bottom 

The simplest model to begin with in discussing underwater gravity correc- 
tions is that of a flat, horizontal bottom at depth h below the sea surface 
and a difference between the sea surface and the reference spheroid (approxi- 
mated by mean sea level) of a vertical distance (h - e) (a tidal correction) 
(Fig. 1). In this model no terrain correction is necessary due to the absence 
of bathymetric relief. 

The first step in correcting the observed gravity is to remove the gravi- 
tational effect of the water above the meter. This is accomplished by calcu- 
lating the effect of a ’Bouguer plate’ of water of density and thickness h. 
The general equation for this correction (Ag^^) is: 

Ag- = 27ra Gh, (1) 

®1 w 

—8 

where G is the universal gravitational constant (6.670 x 10 cgs units). 



4 



For a water density of 1.027 g/cm 
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and h measured in meters, 



eq. 1 reduces to: 



Ag^ = (0.0430 h) mgal. 



( 2 ) 



Since the water is attracting the gravimeter beam or weight upward, this cor- 
rection is added to the observed gravity. 

The second step in the correction requires, in essence, bringing the meter 
from the ocean bottom to the reference spheroid using the free-air correction. 
The general equation for this correction (^§ 2 ^ * 

Ag 2 = 2 GMe/R^, (3) 

where e is the depth of the bottom below the reference spheroid (mean sea 
level) (Fig. 1), M is the mass of the Earth (5.976 x 10^' g) , and R is the 
Earth's radius (6371 km) (MacDonald, 1966), or: 

Ag^ = (0.3083 e) mgal, (4) 

where e is measured in meters. This correction is subtracted from the observed 
gravity as the meter is always below the reference spheroid. The value for 
(h - e) is numerically positive for tidal heights above mean sea level and 
negative for heights below this level. 

The third and final step requires filling in the ocean (now air between 
the reference spheroid and the bottom) with rock if these gravity observations 
are to be correlated with land measurements. This step is numerically analo- 
gous to the 'Bouguer' correction for sea-surface gravity surveys. The equa- 
tion for this correction (Ag^) is: 

Ag^ = Zira^Ge, (5) 

where a , the rock density for the Bouguer plate, is often taken to be 2.67 
g/cm^ (Bobrin, 1960). Grant and West (1965, p. 240) provide a discussion on 
the proper selection of the value for the Bouguer plate density as it is impor- 
tant to relate this value to the local geology in order to obtain a good corre- 
lation with land gravity. Using a rock density of 2.67 g/cm^ (a value that will 
be used throughout this paper), eq. 5 reduces to: 

Ag 3 = (0.1119 e) mgal, (6) 

where e is measured in meters. This correction is added to the observed gra- 
vity as mass is being added below the reference spheroid. 

Combining all these corrections into one expression (Ag. ) gives: 

b 
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(7) 



= Agj^ - Ag^ + Ag^ 



= 2 t[o Gh + 2G(7Ta 
w r 



- M/R^)e, 



or : 



Ag^ = (0.0430h - 0.1964e) mgal. (8) 

Since the tidal height (h - e) will be generally small relative to h, this total 
correction will usually be negative. 

Because of the absence of bottom relief for this model, eliminating the 
necessity for a terrain correction, the Bouguer anomaly (BA) is defined by the 
relationship : 

BA = g^ + Ag^ - g^, (9) 

where g is the gravity value observed by the underwater meter and corrected 
for ear?h tides, instrument drift, and earth curvature, and tied to a reference 
point within the international network of absolute gravity values. The term g , 
the theoretical value for gravity on the surface of the reference spheroid, is^ 
given as a function of latitude (X) by the International Gravity Formula 
(Heiskanen and Vening Heinesz, 1958): 

g^ = [978.0490 (1 + 0.0052884 sin^X - 0.0000059 sin^2X)] gal. (10) 

Bouguer anomalies as calculated from eq. 9 for underwater gravity stations may 
be correlated with Bouguer anomalies for land stations. 

Sea-surface gravity measurements, when corrected for horizontal and vertical 
ship accelerations, cross-coupling effects due to interaction between horizontal 
and vertical accelerations on the meter, instrument drift, earth and ocean tides, 
earth curvature, and the Eotvos effect, give free--air gravity anomalies when com- 
pared with g as formulated in eq. 10. In order to correlate these anomalies 
with underwater gravity measurements, a slightly different set of corrections to 
the observed bottom value is necessary. The first step is to remove the upward 
attraction of the overlying water using Ag from eq. 1. Next, the meter reading 
is reduced from the bottom value to the value on the reference spheroid through 
use of Ag 2 from eq. 3. Finally, the ocean is again filled with water below the 
reference spheroid using the following correction (Ag^) : 

Ag, = 2710 Ge, (11) 

4 w 

or : 

Ag^ = (0.0430 e) mgal, (12) 
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where e is measured in meters. This correction is positive as the mass of 
water added increases the downward gravitational attraction on the meter beam 
when located on the reference spheroid. The combination of Ag and Ag is called 
the double Bouguer correction (Sheriff, 1973, p. 61). Combining these^^correc- 
tions into one expression (Ag^) gives the relationship: 

Agf = Ag^ - Ag^ + Ag^ 

= 27Ta^Gh + 2G(iTa^ - M/R^)e. ( 13 ) 

Using values for a^, G, M, and R as defined previously and measuring h and e 
in meters, this equation reduces to: 

Ag^ = (0.0430h - 0.2652e) mgal. (14) 

As with the expression for Ag in eq. 8, this correction will usually be negative. 
Thus, the free-air anomaly (FaA) for the underwater gravity measurement correla- 
table with sea-surface free-air anomalies is defined as: 

FAA = + Ag^ - g^. (15) 



Smooth sloping bottom 

In this model (Fig. 2), the terrain correction is introduced to reduce the 
underwater gravity measurements. Although the bottom is sloping, there are no 
bathymetric or topographic irregularities in the model. The most logical approach 
to this problem is first to make terrain corrections that result in a model for 
a flat, horizontal bottom, and then to reduce the gravity data in the same manner 
as was discussed in the previous section. Because the bottom is immersed in 
water of finite density rather than in air with a negligible density and because 
some of the ’bottom’ may be above the reference spheroid due to the close proximity 
of the shoreline, certain modifications of the weighting factors for the zones 
in the terrain correction graticule are necessary. The terrain correction chart 
and tables developed by Hammer (1939) may be used" these tables assume a density 
of 2.0 g/cm^ for the material making up the relief. Use of digital computers 
in place of the charts and tables can significantly reduce the time required to 
make these corrections (Kane, 1962; Bott, 1959). 

The first step in applying the terrain correction to the underwater station 
in Fig. 2 is to fill in segment A (a wedge of water between the station depth 
and the bottom seaward of the station) with rock of density a . Since water of 
density a already fills this wedge, the correction density of the filled wedge 
has the value (a - a ). Assuming densities of 1.027 g/cm^ and 2.67 g/cm^ for 
a and a , respectiveYy , and using Hammer’s correction tables, the total terrain 
correction for compartments in this segment is multiplied by a factor of 0.822 
[i.e. , (a - a )/a , where a is the density used to make up the selected 
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correction table] . This correction is added to the observed gravity as it adds 
mass below the meter depth. 

The next step requires removing the effects of the wedge of rock above the 
station depth (segments B and C in Fig. 2). There are two parts to this step 
because the material in segment B is immersed in water, while in segment C it is 
immersed in air. The total terrain correction for segments B and C is first 
calculated assuming a density of a^; again using Hammer’s tables and assuming a 
value of 2.67 g/cm for a , this correction must be multiplied by a factor of 
1.335 (i.e. , • This part of the correction is also added to the observed 

gravity as it represents removal of mass above the meter location. At this point 
it is obvious that too much mass has been removed by this correction because in 
segment B part of the terrain effect of the rock is balanced by water of finite 
density. Thus, the second part of this step of the terrain correction involves 
calculating the effect of filling in segment B with water. The total terrain 
correction for this segment must be multiplied by a factor of 0.514 (i.e., 
a /a ) if using Hammer’s tables. This part of the correction is subtracted from 
t!ie observed gravity as it represents the addition of a slab of water of thickness 
h (segment B) above the meter depth. 

In summary, the terrain correction involves filling in segment A with rock 
rather than water, filling in segment B with water rather than rock, and re- 
moving the effect of the rock above the sea surface in segment C. For the simple 
model demonstrated by this case, Sandberg (1958) gave tables of terrain correc- 
tions for an inclined plane with slope angles of 1^ to 30® based on Hammer’s 
graticule zones and density. If the actual bathymetry and topography of the 
study area can be approximated by a uniform slope, these tables may be used as 
outlined above for Hammer’s tables. 

The model for this case is now exactly the same as for the previous case, a ' 
flat, horizontal bottom, and the gravity data may be corrected to a free-air 
(eq. 15) or Bouguer (eq. 9) anomaly as outlined previously. Grant and West 
(1965, p. 240) state that ’’unless the heights (or depths) of the topographic 
features exceed one-twentieth of their distance from the reference point, they 
are usually omitted”. This would mean for this case that the terrain correc- 
tion is probably unnecessary for a bottom slope of 3® or less. 

To correlate the gravity values taken over the continental shelf from the 
underwater meter with sea-surface gravity measurements, it is probably best to 
first make terrain corrections to the bottom data before calculating the free- 
air anomaly using eq. 15. This statement is based on the concept of downward 
continuation (Peters, 1949; Trejo, 1954) which shows that for a relatively 
horizontal bottom the terrain effect on a bottom gravimeter is greater than on 
a meter at the sea surface at some distance above the same station. In shallow 
water where the vertical angle from a given portion of bath 3 rmetric relief to 
the bottom station may be greater than the similar angle from the feature to 
the surface meter, the terrain effect may be greater for the surface meter, 
requiring terrain corrections to be made to the data from both meters. 
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Irregular bottom 



Although this case illustrates the usual model found on the continental 
terrace, the previously discussed cases are applicable approximations in many 
areas. Along the central California coast in the vicinity of Monterey Bay the 
presence of submarine canyons which begin just seaward of the shoreline and 
the rugged coastal mountain range cause a rather narrow continental shelf (about 
10 km wide) and a very steep continental slope (in places, bottom slopes of up 
to 45° have been found) . Rock outcrops and fault scarps on the shelf cause 
many bathymetric irregularities. The model for this case may be broken into 
the same three segments that were discussed in the previous case, and the general 
method for making terrain corrections is the same, although Sandberg’s correc- 
tion tables for an inclined plane are probably not appropriate. Just south of 
Monterey Bay, Souto (1973) calculated terrain corrections of 5.5 to 5.7 mgal 
for stations at water depths of 60 to 75 m near the Carmel Submarine Canyon, 
a tributary of the Monterey Submarine Canyon. Inclinations of 28° were found 
on the continental slope in this area. With Bouguer anomalies ranging from 
+23 to +30 mgal for this area, the terrain correction is an important factor 
in data reduction. 



CONCLUSIONS 

Correcting underwater gravity data presents some unique problems to an 
investigator. If terrain corrections are necessary, it appears to this writer 
that they should be applied first to the observed gravity data to reduce the 
bath 3 rmetry to a model of a flat, horizontal bottom. Next the upward attrac- 
tion of the water above the station should be removed and the station should 
be relocated on the reference spheroid with a free-air correction. Finally, 
water should be again added to the ocean beneath the spheroid if a free-air 
anomaly chart is desired for correlation with sea-surface gravity data, or 
rock should be added if a Bouguer anomaly chart is required for correlation 
with land gravity surveys. This order of application of corrections should 
make more apparent the physical meaning of each step and aid in determining 
which steps, if any, may be simplified by approximation for special bathymetric 
cases or eliminated. 
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Fig. I. Model for Flat, Horizontal Bottom. 



10 




Fig. 2. Model for Smooth Sloping Bottom with Segments for Terrain Correction 
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